Draft version January 13, 2013 

Preprint typeset using I^'T^]X style emulatcapj v. 11/10/09 



THE HIGHLY DYNAMIC BEHAVIOR OF THE INNERMOST DUST AND GAS IN THE TRANSITION DISK 

VARIABLE LRLL 31 

Flaherty, K.M. \ Muzerolle, J. ^, Rieke, G. \ Gutermuth, R. Balog, Z^, Herbst, W. ^, Megeath, S.T. ^, 

KuN, M.^ 

Draft version January 13, 2013 



o 



< 

o 

6 



00 

p 

rn 
O 



X 



ABSTRACT 

We describe extensive synoptic multi- wavelength observations of the transition disk LRLL 31 in the 
young cluster IC 348. We combined four epochs of IRS spectra, nine epochs of MIPS photometry, 
seven epochs of cold-mission IRAC photometry and 36 epochs of warm-mission IRAC photometry 
along with multi-epoch near-infrared spectra, optical spectra and polarimetry to explore the nature 
of the rapid variability of this object. We find that the inner disk, as traced by the 2-5/im excess 
stays at the dust sublimation radius while the strength of the excess changes by a factor of 8 on 
weekly timescales, and the 3.6 and 4.5/im photometry show a drop of 0.35 magnitudes in one week 
followed by a slow 0.5 magnitude increase over the next three weeks. The accretion rate, as measured 
by Pa/3 and Br7 emission lines, varies by a factor of five with evidence for a correlation between 
the accretion rate and the infrared excess. While the gas and dust in the inner disk are fluctuating 
the central star stays relatively static. Our observations allow us to put constraints on the physical 
mechanism responsible for the variability. The variable accretion, and wind, are unlikely to be causes 
of the variability, but are both effects of the same physical process that disturbs the disk. The lack 
of periodicity in our infrared monitoring indicates that it is unlikely that there is a companion within 
~ 0.4 AU that is perturbing the disk. The most likely explanation is either a companion beyond ~ 0.4 
AU or a dynamic interface between the stellar magnetic field and the disk leading to a variable scale 
height and/or warping of the inner disk. 



1. INTRODUCTION 

The widely accepted view is that a pre-main sequence 
star slowly evolves from a deeply embedded core to a fully 
revealed star. After initial collapse of a dense core, the 
star enters the class O/I phase in which it is rapidly ac- 
creting mass from its massive envelope and disk. As the 
envelope dissipates and the accretion rate decreases, the 
star enters the class II phase where it becomes optically 
revealed and is being fed slowly by a viscously-evolving 
irradiated accretion disk. The dust in the disk eventually 
coalesces into large grains and the gas dissipates from the 
system leaving a class HI star. 

Confounding this picture of slow evolution, there is 
growing evidence that these systems are in fact highly 
dynamic. It has been known since the discovery of 
T Tauri stars, which are at the heart of class II sys- 
tems, that the star is usually variable (|Jovlll945[ ). Hot 
and cold spots rotating across the surface of the star, 
as well as variable accretion onto the star can create 
changes of a ma gnitude or more in only a few days in 
the optical fe.g. iHerbst et all I1994D . Recent observa- 
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tions have shown that the circumstellar disk, as traced 
by the infrared emission, can also change dramatically 
in only a few days. T he shape of the disk contin- 
uum appears to ch ange ("Morales-Calder on et al.l 120091 : 
lJuhasz et all 120071 : Muzerolle et al. 2009) and the sil- 
icate feature, which traces the properties of the dust 
in the upper laye rs of the disk c l ose to the star 
also varies ra pidly (|Barv et al.l 120091 : iSitko et al.l 120081 : 
iSkemer et aT1l2010( ). The behavior may be related to 
large structural changes in the inner disk (Sitk o et al.l 
12008"), large accretion events, such as in FU Ori outbursts 
(iHart m ann fc Kenvon. 1996,1 or to turbulence within the 
disk (jTurner et al.ll2010D . The previous studies provide 
some insight into the source of the variability, but are 
limited by the lack of ancillary data and limited tempo- 
ral or wavelength coverage. 

To understand better what could be causing these 
changes, we have studied one object, the highly variable 
star LRLL 31, in detai L This object shows dramatic 
mid-infrared variability (jMuzeroUe et al.ll2009f l that may 
be due to a non-axisymmetric s tructure in the inner 
disk with a changing scale height (jFlahertv fc Muzerolle! 
120101). We report continued monitoring in the mid- 
infrared, as well as multiple epochs of near-infrared spec- 
tra, optical spectra and polarimetry measurements to un- 
derstand the dynamics of the system better. Using these 
observations, we study the changes in the dust, the gas 
and the central star, and any connections among them. 
The near simultaneity of these data sets allow us to rule 
out specific mechanisms that might be responsible for the 
observed variability. 



2. DATA 
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Our monitoring of LRLL 31 consists of multiple epochs 
of ground and space based observations covering the op- 
tical as well as the near and mid-infrared. We obtained 

0- S-Sj um spectra with the SpeX instrument (jRavner et alJ 
[2OOI on the NASA Infrared Telescope Facility (IRTF) 
in fall 2008 and 2009. Table [D lists the observing log 
for all of our observations, and the spectra are shown 
in Figures [1] and [2l On all but two nights we used 
both the SXD (0.8-2. 2/im) and LXD2.1 (2-4.8Mm) grat- 
ings, while on October 8,9 2009 we only obtained SXD 
spectra. On Oct 10,18,19 2008 we used the 0.5"xl5" slit 
(R^1500) while on the other nights we used the 0.8"xl5" 
(R^800) slit due to poorer seeing conditions. We also 
supplem ent our observations with SXD spectra take n 
in 2005 (|Muzerolle et al.|[200 9h and 2006 fDahm"200i). 
The data were reduced using Spextool (Gushing ct al. 
120041) and corrected for telluric absorption using xtell- 
cor (jVacca et al. 2003). Our telluric standards were the 
nearby AOV stars HD 19600 and HD 23441. The to- 
tal exposure times were 30 minutes for the SXD and one 
hour for the LXD spectra. A set of calibration exposures, 
including flat field and argon arc spectra, were obtained 
for each target-standard star pair. 

We also observed the bright weak line T Tauri star 
(WTTS) HD 283572, which shown no sign of accretion or 
an infrared excess, with an exposure time of 10 minutes. 
It has the same spectral type (G6) as LRLL 31 and is 
used as a comparison star when measuring the veiling 
and strength of the emission lines. We took SXD spectra 
of this star with a slit of 0.8"xl5" on 5 Jan 2010. The 
spectra were reduced in the same manner as the LRLL 
31 spectra using the telluric standard HD 25152. 

In addition to the spectra taken by Spex, we also ob- 
tained JHK photometry (Table [2]) from the slit-viewing 
near-infrared guide camera. LRLL 31 was observed with 
three dithers, standard reduction techniques were used 
and aperture photometry was performed on the two stars 
in the field of view. We did not observe any photo- 
metric standards, but instead derived the photometry 
of LRLL 31 relative to the other star in the field (LRLL 
169, 2MASS03441776-I-3204476). This allows us to de- 
rive accurate relative photometry, but there could be 
large uncertainties in the absolute photometry. LRLL 
169 is a diskless cluster member and is not marked 
as variable (g=0.0 6 mag) in the optical monitoring of 
iCohen et all (pOOl . Our Spitzer warm mission photom- 
etry (described below) of LRLL 169 shows variations of 

1- 2% at 3.6 and 4.5/im which is consistent with the un- 
certainty in the photometry. This suggests that the pho- 
tometric variability is no more than a few percent, which 
is smaller than the absolute uncertainty (0.1 mag) asso- 
ciated with using only one standard and we ignore it. 

We also include high-resolution Ha spectra from 
Keck (kindly provided by Scott Dahm) and from the 
MMT. The MMT spectru m was taken with Hectochelle 
(|Szentgvorgvi et al.|[l998[) in 29 Feb 2008 with a resolv- 
ing power of R^ 30, 000. The Keck spectra were taken 
with HIRES C Vogt et all [1991 on 28 Nov 2006 and 03 
Dec 2008 with a resolving power of R~ 33, 000. Ob- 
serving details for the MM T spectra are reported in 
iFlahertv fc MuzeroUd (|2008[ ) while the reduction process 
for the Keck spectra is reported in iDahml ()2008D . 

We include spectro-polarimetric observations 
taken on several nights using the SPOL imag- 



ing/spectropolarimeter () Schmidt et al.l 119921) on the 
Steward Observatory 2.3 m telescope. These low- 
resolution optical spectra, shown in Figure |3l extend 
from 4000- 7500 A, and from them we can extract the 
strength of the Ha line as well as the shape of the 
continuum, which is useful for measuring the reddening 
toward the source. Details on the observa tional mode 
and the data reduction can be found in iSmith et al.l 
(2003). An intermediate resolution optical spectrum of 
LRLL 31 was obtained on 2009 October 11 using the 
CAFOS instrument on the 2.2-m telescope of the Calar 
Alto Observatory. The exposure time was 1880s. The 
R-lOO grism covered the wavelength range 5800-9000 A. 
The resolving power, using a 1.5-arcsec slit, was R^3500 
at 66OOA. We observed the spectrum of a He-Ne-Rb 
lamp for wavelength calibration. The spectrum was 
reduced using standard IRAF routines. 

We have observed LRLL 31 on multiple epochs with all 
three instruments onboa rd the Spitzer Space Telescope. 
IRS (|Hoiick et al.l [200I data were originally obtained 
on 09 Oct 2007, 16 Oct 2007, 24 Feb 2 008 and 02 Mar 
2008 a nd these spectra were presented in lMuzeroUe et aP 
(|2009f ). The IRS data reduction pipeline has been im- 
proved and the updated spectra are presented here. The 
spectra were extracted using the Spectroscopic Mo deling, 
Analy sis and Reduction Tool (SMART; Higdo n et alJ 
(|2004f l) version 8.1.2, starting with the basic calibrated 
data products from the Spitzer Science Center reduction 
pipeline version S18.7. Rogue pixels were first removed 
using the "IRS_CLEAN" program, constructing a pixel 
mask for each of the SL and LL modes using off-order im- 
ages as a reference. SMART was then used with optimal 
manual point source extraction using local sky subtrac- 
tion. The local sky background was calculated by fitting 
a polynomial to all sky pixels on either side of the source 
at each row in the spatial direction. The final calibrated 
spectra for each nod and order were then combined using 
the sigma-clipped averaging function in SMART; typical 
measurement uncertainties were about 2% at most wave- 
lengths. 

IRAC (jFazio et al.ll2004|) observations were taken dur- 
ing the cryo-mission, when all four channels were avail- 
able, on 19-23 March 2009 (PIP 50596) (Ta ble H along 
with the original GT O (|Lada et al.l I200I and C2D 
(| Jorgensen et all l2006( ) data taken on 11 Feb 2004 and 
08 Sep 2004. A warm-mission observing program (PID 
60160) obtained repeated 3.6,4.5/Ltm photometry at a 
varying cadence with observations ranging from every 
four hours to every other day over most of the 40 day 
visibility window of IG 348 (Table fTTI). The IRAQ 
data r eduction pipeline is described in lGutermuth et al.l 
()2009( ) with updates appropriate for the warm Spitzer 
missio n described in G utermuth et al. (in prep). 
MIPS (iRieke et al."2004D observations were taken on 21 
Feb 2 004lLada ct al. 2Q0l), 19 Sep 2004 (jRebull et al.l 
[2007) and 24-27 Sep 2007, 12 Mar 2008, 19 Mar 2008 
(iMuzcroll e et al.ll2009f). Details on the re duction of the 
MIPS data are in iMuzeroUe et al.l ()2009[ ). Information 
on all of these observations is included in Table [T] 

3. EVIDENCE FOR FLUCTUATING GAS AND DUST 
PROPERTIES WITH STATIC STELLAR PROPERTIES 

3.1. Stellar Properties 
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The stellar luminosity can be derived from the J band 
assuming it is dominated by stellar flux, as demonstrated 
in Figure [1] once it has been dereddened. We use the 
shapes of both the near-infrared and spectro-polarimetry 
spectra to measure the reddening of LRLL 31. The spec- 
tra can be dereddened until their shapes match that of 
a standard of the same spectral type. We ignore the 
K band spectrum since it could be modified by emis- 
sion from dust, rather than jus t extinction. We use 
the Ry=5.5 extinction law from iWeingartner fc: Drain^ 
(j2001) in both the optical and near-infrared, which is 
more appropriate for a star embedded in a molecular 
cloud than Ry=3.1, and provides a good fit in com- 
paring LRLL 31 and the standard. For our standard 
we use the SpeX spectrum of HD 283572 in the near- 
infrared, and in the op tical we used a Kurucz model 
(Teff=5750, log g=2.5) (jCastelli fc Kuruczl [2001 . The 
results are listed in Table |31 Since the shapes of the 
standard and the spectral do not match perfectly due to 
effects such as imperfect telluric correction and scattered 
light within the instrument, there is a 0.5 mag uncer- 
tainty in the best fit extinction. We find that on average 
the extinction is Ay ~ 8.8. This i s lower than previ- 
ous estimates ([Muzerolle et al.l 12009^ because of our use 
of a different extinction law, which will convert the same 
amount of reddening to a smaller value of Ay. Occasion- 
ally the extinction appears to drop down to Ay ^ 8.2, 
but these values are consistent with the average given 
the large error bars. Figures [T] and [5] show the observed 
near-infrared spectra compared with a reddened stan- 
dard along with the observed JHK photometry and si- 
multaneous 3.6,4.5/im photometry where available. The 
agreement between the shape of the spectra and the pho- 
tometry suggests that we are accurately measuring the 
shape of the SED. Figure [3] shows the optical spectra 
along with a reddened standard. Using the optical spec- 
tra we derive extinctions that are consistent with the 
near-infrared, including occasional apparent drops in Ay. 

The reddening is combined with the J band photom- 
etry, a bolometric correction from Kcnvon & Hartmann 
(11995!) and a distance of 320 pc (|Luhman et al.ll2003l) to 
produce the stellar luminosity, which is listed in Table O 
The uncertainties are due to the the errors in the photom- 
etry and the extinction and are typically O.GLq. There 
is a range of luminosities, but given the error bars they 
are consistent with a constant luminosity of A.SLq. This 
luminosity corresponds to a radius of 2.1Rq given a G6 
(Te//=5700) spectral type. If the low values of luminos- 
ity represent a real deviation from the mean, they may 
illustrate moments where the star is partially obscured 
by dust that is not taken into account in our estimate 
of the extinction. Overall the star appears to be un- 
derluminous relative to a G6 star at the age of IC 348 
(2-3 Myr), based on the lSiess et al.l (|2000( ) isochrones, al- 
th ough it is roughly consi stent wi th the 3 Myr isochrone 
of IPalla fc StahleJ (fl999l) which ILuhman et all (l200l 
finds is more appropri ate for the high m ass members of 
IC 348. Based on the ISiess et all (|2000D isochrones and 
the observed luminosity, LRLL 31 is roughly 1.5 Mq. 

Recent monitoring in I-band (Baliber et al. in prep) 
shows small variations of ~ 3% with a period of 3.4 
days, which are consistent with the rotation of cool 
spots across the stellar surface suggesting that 3.4 days 
is the rotation period of the star. Despite the fre- 



quent monitoring of I C 348, thi s periodicity has not 
been observed (Cohen erall 120041: iLittlefair et "all 120051: 
iCieza fc Baliben ,2006,: iNordhagen et al.l I2006D . These 
previous observations were ta ken in 1998-2004 and found 
1=15.24 C Cohen et all 120041 ) and 1=15.36 with occa- 
si qnal dips down to 1= 15.5 (see star 0130 in Figure 4 
of ILittlefair et all 1200 5V while the most recent observa- 
tions measured 1=15.82 for LRLL 31. The near-infrared 
photometry shows a long term decrease from the 2MASS 
data (operated from 1997-2001) to our measurements in 
2008 and 2009. This suggests that there may be a long 
term decrease in the stellar fiux. Continued optical mon- 
itoring of this star will be useful in confirming the peri- 
odicity and dimming of LRLL 31. 

Our J-band photometry does not show a large change 
in the stellar flux, and some of the optical monitoring 
by Baliber et al. overlaps the first week of our Spitzer 
warm-mission monitoring and does not show a substan- 
tial decrease in the optical fiux, which leads us to as- 
sume that the stellar flux is constant during our infrared 
monitoring. With the per i od in formation and the vsini 
measurement from lDahml (|2008l ) we can estimate the in- 
clination of this star to be i = 38 ± 7°, which includes 
both the uncertainty in the v sin i measurement and the 
radius as derived from the luminosity. This is the incli- 
nation to the stellar rotation axis, and it is possible for 
the stellar magnetosphere and the disk to be misaligned 
with this axis. 

We find a change in the radial velocity based on three 
epochs of high resolution optical spectra. On 28 Nov 
2006 the radial velocity was 12.82±0.42 km/sec, on 29 
Feb 2008 it was 8.36±0.27 km/sec and on 03 Dec 2008 it 
was 15.7±1.0 km/sec. These changes could be due to the 
infiuence of a companion, or partial obscuration of the 
surface by cool spots or circumstellar material. If part of 
the stellar surface, such as the side rotating toward us, 
is partially blocked then the radial velocity will appear 
to become red-shifted. As the obscuring material rotates 
around the star the apparent radial velocity will change. 
The rotation of a spot across the surface of the star has 
been sugg ested to explain the r adial velocity variations in 
AA Ta u (.Bouvier et aT]l2007af) and the transitional disk 
T Cha (|Schisano et al.ll2009() . The rotational velo city of 
LRLL 31 has been measured to be 19.8 km/sec (|Dahml 
|2008() . which is consistent with the observed radial ve- 
locity variation, although the size of the spot required 
to produce such large radial velocity variations is bigger 
than typically seen around G-type stars. A change in 
the shape of the bisector of the cross correlation func- 
tion used to derive the radial velocity would indicate the 
presence of obscuring source rotating in front of the star, 
but we do not have the signal to noise to test this theory. 
A clump of obscuring material, either in the disk or in 
the molecular cloud, would need to block almost half of 
the star in order to create the observed radial velocity 
variation. This would produce a drop in the optical flux 
of 0.75 mag, along with a substantial increase in the ex- 
tinction, neither of which are observed. If these radial 
velocity measurement are confirmed then the most likely 
explanation is the presence of a companion somewhere in 
the system, although this does not prove that the com- 
panion is responsible for the observed infrared variabil- 
ity. In section 4.3 we place constraints on the mass and 
location of such a companion using the available data. 
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although more radial velocity measurements are needed 
to confirm and fully characterize this companion. 

3.2. Dust Properties 

LRLL 31 was originally singled out based on its mid- 
infrared variability (Muzerolle et al-- 2009). In IRS spec- 
tra separated by one week the 5-8/im flux dropped al- 
most to the photosphere while the 8-40^m flux increased. 
This behavior was unique and difficult to explain given 
typical sources of variability, such as the rotation of 
hot/cold spots across the surface of the star or extinc- 
tion events. Modeling finds that the wavelength depen- 
dence, and timescale, of the variability is consistent with 
fluctuations in the scale height of the inner disk on a 
dynamical time scales, which is on the order of one week 
([Flaherty fc M uzerolle 2010). As the inner disk grows 
it increases the short wavelength fiux, while at the same 
time shadowing the outer disk, decreasing the long wave- 
length flux. The previous observations of the continuum 
emission were limited to probing material far from the 
star, at 1-10 AU, which may not be the source of the 
variability. The silicate emission may arise from closer 
to the star (e.g. [Kesslcr-Silacci et al. 2007), but this lo- 
cation is difficult to constrain without detailed radiative 
transfer models, especially given the transition disk na- 
ture of LRLL 31 

Our 3-5/im spectra give us information on the dust 
close to the star. By subtracting the stellar photosphere 
from the spectra we are left with only the emission from 
the dust. When the dust at the inner edge of the disk 
is optically thick, the emission should appear as a single 
tempe rature blackbody at the dust destruction temper- 
ature ([Muzerolle et al.ll2003l: lEspaillat et al.ll2010[) . The 
shape of the excess spectra will tell us the temperature of 
the hot dust which translates to the location of the inner 
edge of the dusty disk. To subtract the photosphere we 
must first normalize the standard to the LRLL 31 stellar 
component. We use the veiling at 2.15/im to separate 
the stellar and excess emission of LRLL 31. 

We can derive the veiling based on the strength of pho- 
tospheric lines in the K band relative to the G6 standard. 
Vefling is defined as r = F/Fphot where F is an excess 
continuum flux that comes from the accretion disk; as r 
increases the strength of the lines relative to the contin- 
uum decreases. The value of r can be used to measure 
the amount of excess emission. We calibrate the veiling 
by adding a constant flux to our G6 standard in 16 small 
bands 0.05/im wide from 0.8-2.5/im and comparing the 
veiled standard to the observed spectra. The best-fit is 
determined independently in each band using the mini- 
mum chi-squared. These multiple measurements provide 
r as a function of wavelength and we use this linear fit to 
derive the veiling at 2.15/im. There is substantial scat- 
ter between individual measurements, but we are able to 
use the line fit to estimate r with an uncertainty of ±0.1. 
Figure 131 shows the measurements in our 16 bands for one 
day, demonstrating that we can obtain good fits across 
much of the spectrum, showing a trend of increasing veil- 
ing with increasing wavelength. The results are listed in 
Table ISl 

Once we have used the veiling to determine the excess 
emission at 2.15/xm we normalize the standard and LRLL 
31 spectra at this wavelength. Since our G6 WTTS stan- 
dard spectrum only extends to 2.3/Ltm we use a Kurucz 



model (Teff=5750, log g =2.5) ([Castelh fc Kuruczll200l 
to extend it out to 5/Ltm. The standard and LRLL 31 
are subtracted yielding an excess spectrum whose flux is 
relative to the photospheric flux at 2.15/im. These excess 
spectra are shown in Figure (5] 

To find the temperature of the hot dust we fit a sin- 
gle temperature blackbody to the excess spectrum. We 
artificially increased the uncertainties in the wavelength 
ranges 2.4 < A < 3.5/im and A > 4.3/im in order to ac- 
count for systematic effects from an incomplete telluric 
correction, which affect some regions of the spectrum 
more than others. This is especially important for 08 
Nov 2009 when the large airmass difference between the 
telluric standard and LRLL 31 (0.48) resulted in a poor 
correction of the telluric features. Typical uncertainties 
on the derived temperature are 250K, which is mostly 
due to the uncertainties in the veiling and in the extinc- 
tion. As seen in Figure [5] a blackbody is a good fit to 
the excess spectra. The full spectra, along with the veil- 
ing measurements as a function of wavelength are shown 
in Figure [6] 

We can double check the shape of the excess spectra by 
normalizing both the standard and LRLL 31 at 1.15/im, 
where we expect and measure zero veiling. This avoids 
using the veiling to normalize the spectra, since the large 
uncertainties in the veiling lead to large changes in the 
derived temperature. Veiling presents an extinction in- 
dependent measure of the shape of the LRLL 31 spectra 
relative to the standard, which can be used in this case to 
confirm the shape of the spectra, rather than relying on 
them to define the shape. We find that the temperatures 
derived in this instance are consistent with our previous 
derivation within the uncertainties. The largest discrep- 
ancies are for the 31 Oct 2009, where the new tempera- 
ture is 1900 K, and 11 Oct 2008, where the new tempera- 
ture is 1620 K. Normalizing the spectra at 1.15/(m brings 
these two measurements into better agreement with the 
derived temperatures from the other five nights. 

The temperatures, listed in Table|6l are consistent with 
each other, within the uncertainties, and the average 
temperature of 1830 K is consistent with the dust sub- 
limation temperature (^e.g. iDullemond fc Monnier 2010) . 
This implies that the sublimation of dust sets the inner 
edge of the disk rather than some other process, such as 
the sculpting of the disk by a companion. Converting the 
dust temperature to a location in the disk depends on the 
properties of the dust, the luminosity of th e star, and also 
on the covering fraction of the inner rim ([Isella fc Nattal 
l2005l) . If one of these dependencies was rapidly varying 
then the location of the inner rim could fluctuate while 
the temperature stays constant. Our measurements indi- 
cate that the accretion luminosity is much less than the 
stellar luminosity, the stellar luminosity docs not vary 
significantly, and we do not expect the dust properties, 
such as grain size and crystallinity, to change on the ob- 
served timescales, suggesting that a constant dust tem- 
perature corresponds to a constant sublimation radius 
(- 0.15AU). 

While the average temperature of the hottest dust 
is roughly consistent with the dust sublimation tem- 

^ We derive the extinction assuming that there is very little 
excess in the J and H bands, which could potentially change the 
shape of the spectra. 
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perature, it is s till hi gher than is typicall y measured 
(|Muzerolle et alj 120031: [Espaillat et aI1l2010D . This dis- 
crepancy may either be due to unique gas and dust prop- 
erties that allow for higher than normal temperatures 
or systemic effects in our measurements that cause us 
to overestimate the dust temperature. Dust composed 
of irons and olivines at high gas densities has a higher 
sublimation temperature than dust particles more com- 
monly seen around young stars (jPollack et al.llT99^ . An 
inner disk composed mostly of this material would have 
an abnormally high temperature. While we cannot di- 
rectly measure the dust and gas properties of the inner 
disk we can use the silicate feature and the accretion 
rate as proxies. The silicate feature is sensitive to dust 
properties at radii of a few te nths of an AU to a few 
AU (Kcssl er^Silacci et al.ll2007[ ) and the dust properties 
responsible for this feature generally do not change by or- 
ders o f magnitude throughout the disk ( van Bockcl et all 
12004) . The shape of the silicate feature in LRLL 31 is 
consistent with typical ISM-like silicate dust with sub- 
micron sized grains, not a heavy mix of irons and olivines. 
Also the accretion rate, which is proportional to the gas 
density in a viscously evolving disk, is inconsistent with 
an extremely high gas density. The fact that this star is 
a transition disk would suggest that there is a decrease 
in the gas density rather than an increase. 

This suggests that there is a systematic effect in our 
measurements that lead to an overestimation of the dust 
temperature. Overestimating the veiling or the extinc- 
tion, possibly due to an imperfect match between the 
standard star and the intrinsic photospheric emission of 
LRLL 31 will lead to a higher measured dust tempera- 
ture. Given the difficulty in measuring the veiling at such 
low levels, this is the most likely source of the discrep- 
ancy. Our conclusions about the change in veiling should 
be unaffected by this systematic uncertainty since we use 
the same standard, accounting for changes in resolution, 
and derive the veiling using the same method on each 
night. 

While the temperature of the dust is constant, the 
strength of the excess fluctuates throughout our observa- 
tions. The variations in the strength of the excess occur 
on weekly to monthly timescales, but do not show daily 
fluctuations (Figure ^ . We integrate a blackbody at the 
derived temperature and strength for each night to de- 
termine the luminosity of the excess relative to the stellar 
luminosity, which is an estimate the covering fraction of 
the inner disk, since the flux from the disk is reradiated 
stellar flux. Our results are listed in Table H) For the 
nights of 8,9 Oct 2009 we assume the dust is at 1830 K, 
the average of the other seven nights, and fit the black- 
body to the IRAC data to determine the strength of the 
excess. We find that the covering fraction varies by a 
factor of 8, from 1-8%. 

In a typical T Tauri sta r, with a puffed inner r im, the 
covering fraction is 12% (jPullemond et al.ll2001[ ). while 
typical debris disks h ave covering fractions of <0.1% (e.g. 
iPlavchan et"ani2009| ). On average the covering fraction 
of the dust in LRLL 31 is below that expected for an 
optically thick puffed inner disk. In section 4 we discuss 
physical models that can produce changes in the infrared 
excess on short timescales, but here we consider what the 
small average covering fraction implies about the inner 
disk structure. The covering fraction for an optically 



thick disk is proportional to the ratio of disk height to 
radius, H/R (jDullemond fc Monnieii[20lol ). and if we as- 
sume that R stays constant, then we are observing a 
reduction in the scale height relative to a normal disk. 
This deviation in the average scale height from a typical 
disk could be due to (1) a dramatic decrease in the lo- 
cal surface density or (2) a decrease in the dust opacity, 
possibly due to the growth of the dust grains. Both of 
these effects would reduce the r = 1 surface of the disk 
where stellar photons are absorbed, reducing the cover- 
ing fraction. The accretion rate, which in a viscous disk 
is proportional to the surface density, is not unusually 
low for a T Tauri star, suggesting that enhanced grain 
growth and settling may be causing the decrease in the 
covering fraction, rather than a decrease in the surface 
density. 

This long term evolution of the grain size and settling 
toward the midplane may or may not be related to the 
rapid fluctuations seen in our data. The change in cov- 
ering fraction implies a change in the scale height of a 
factor of 8 if the perturbation is axisymmetric. A model 
of the infrared fluctuations does not need to produce ax- 
isymmetric perturbations of the disk but does need to 
create features that are large enough to match the ob- 
servations. If the change in covering fraction is due to a 
localized blob being launched from the disk and the disk 
always has a covering fraction of 1% then the other 7% 
can be due to a blob with a circular diameter of ~50° as 
seen from the star. 

The 3.6 and 4.5/im flux (Figure [7]), where we have more 
complete monitoring, shows a slow variation over the 
course of weeks, consistent with the change in the in- 
frared excess measured from the spectra. The infrared 
flux is a mix of stellar and disk emission and our excess 
spectra taken during the warm mission monitoring sug- 
gest that most of the change in the infrared flux comes 
about from a change in the disk emission, not from vari- 
ability in the stellar emission. Some of the optical moni- 
toring of Baliber et al. overlaps with our infrared moni- 
toring, and the stellar flux does not follow the same trend 
as the infrared data, supporting our interpretation that 
the dust is varying while the star is constant. 

In our mid-infrared observations the flux varies contin- 
uously (Figure [5]). This suggests that the physical cause 
is not a single event, but a continuously occurring phe- 
nomenon. Unlike a FU Ori or EXor outburst, the vari- 
ability in LRLL 31 does not occur once and slowly de- 
cay. Our IRS spectra separated by six months show that 
the wavelength dependence of the variability does not 
change. This suggests that the same mechanism domi- 
nates the variability over many years. It also supports 
our assumption that the variability we study in any one 
year has the same characteristics as in any other year. 

The extinction and polarization (Tables H] and [7|) can 
provide clues to the structure of the disk, assuming that 
they come from the disk. If the grains in the molecu- 
lar cloud are aligned, they can create polarized flux from 
the star. Starlight that is parallel to the long axis of 
the grains will be more highly extincted than starlight 
that is parallel to the short axis of the grains. This 
differential extinction will lead to polarized light leav- 
ing the cloud, with the polarization being proportion 
to the extinction. Two other nearby cluster members 
were observed by SPOL and have previous polarimetry 
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measurements. LRLL 169 (11" from LRLL 31) has a 
P=4.2% and q=6.9° while LRLL 4 (3.5' from LRLL 31) 
has P=2.8% and q=16 6.5°. and the extincti on is Ay=2.8 
and 1.8 respectively (jLuhman et al.l [20031 ). If the po- 
larization is from the molecular cloud then the ratio of 
polarizatio n to extinction should slightly decrease with 
extinction (jGerakines et el.lll995[) . For LRLL 169,4 and 
31 the ratio is 1.5, 1.6 and 0.9 %/mag respectively, which 
is consistent with the expected decrease. The molecular 
cloud density needs to be highly variable in order to ex- 
plain the observations since LRLL 169 is only 3500 AU 
away from LRLL 31 and shows a much smaller extinc- 
tion and a different polarization angle. The high polar- 
ization/extinction could also be coming from the outer 
disk, if the disk is close to edge on and the system is mod- 
erately flared. Differential extinction can work with ma- 
terial in the outer disk in much the same way as it would 
operate in the molecular cloud. An edge-on disk is consis- 
tent w ith the inclination used by iFlahertv fc Muzerolld 
(|2010D to fit the LRLL 31 SED with models of a warped 
disk. We can exclude the inner disk as the source of the 
extinction and polarization because of the lack of strong 
variability in these quantities, whereas we have observed 
continuous and large fluctuations of the inner disk. 



3.3. Gas Properties 

Next we consider the activity of the gas, which we trace 
through a combination of emission lines from the infrared 
and optical. Our observations of the Pa/3 and Br7 lines 
in the infrare d let us measure the ac cretion rate onto the 
star directly (jMuzerolle et al.l [19981 ). To derive the line 
fluxes, we first subtract the photospheric absorption as 
measured by the G6 WTTS standard. In this case we find 
that using a WTTS, rather than a dwarf or giant spec- 
trum, does better at subtracting the photospheric ab- 
sorption especially on the wings of the line. We normal- 
ize the absorption-subtracted spectrum to the continuum 
around the line. These spectra for each day are shown in 
Figure ini The equivalent widths of the lines, as measured 
from the continuum normalized spectra (Figure |9|), are 
listed in Table IS] Veiling will change the strength of a 
line, especially for Br7 where the veiling is significant, 
and in Table [8] we include the veiling corrected equiva- 
lent width (EW), which is derived by artificially veiling 
HD 283572 and subtracting this new standard from the 
observed spectrum. Using our nearly simultaneous mea- 
surements of J and K band photometry to estimate the 
continuum level, we derive the line fiux. Since we use the 
observed K band photometry, which is a mix of stellar 
emission and veiling flux, we had to account for the veil- 
ing when measuring the EW, resulting in the corrected 
EW being lower than the uncorrected EW. The smaller 
EW, combined with the stellar plus veiling flux of the 
K band continuum produces an accurate estimate of the 
line flux. While ot hers have found sign ificant veiling near 
1/im in some stars (iFischer et al.l[2008[ ). we measure little 
veiling in the J band, and assume it is negligible when 
deriving the flux of the Pa/3 line. Assuming a distance of 
320pc to the star, we can then derive line luminosities. 
These line luminositi es are directly correla ted with the 
accretion luminosity ([Muzerolle et al.lll998D and we can 
derive the accretion rate using the formula: 
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The measured line luminosities and accretion rates are 
listed in Table |9l We cannot derive accurate accretion 
rates from the 2005 and 2006 spectra because we do not 
have simultaneous photometry, but the strong emission 
lines suggest that they were on the high end of the range 
of measured accretion rates. There is a large discrep- 
ancy between the accretion rate derived by Pa/? and Br7 
that is unexpected given that both measure the accretion 
rate onto the star. The uncertainties listed in Table [3] 
only include the uncertainty in the photometry and the 
equivalent width measurements and do not include the 
uncertainty in the conversion between line flux and ac- 
cretion luminosity, or the uncertainties in the mass and 
radius of LRLL 31. The conversion uncertainty is closer 
to a factor of 2, much larger than the typical 10% uncer- 
tainty, and can account for the discrepancy between the 
accretion rate derived by Pa/3 versus Br7. Unresolved 
absorption in the line may also explain some of the ob- 
served differences. Further conclusions drawn about the 
accretion rate are based on the Pa/3 because it is the 
stronger and is less susceptible to veiling, which makes 
it a more reliable tracer of the accretion rate. 

There are large variations in the Pa/3 and Br7 emission 
throughout our observations. The EW ranges from al- 
most zero up to 4.5A. The accretion rates vary by a factor 
of 5. The infrared emission lines demonstrate that the 
accretion is changing constantly with the largest changes 
occurring on weekly, rather than daily, timescales. Ac- 
cretion variability is common among T Tauri stars, and 
the daily level of variability seen i n LRLL 31 is similar 
to that seen in many other stars fflartigan et al.l [19911 : 
IGuUbring et"all [l996[ ) . The change seen in 2009, when 
the Pa/3 line goes from absorption to emission is less com- 
mon. 

Our measurements of the Ha equivalent width {Ta- 
ble [TU|) from the optical spectra show emission at every 
epoch, implying that accretion is continuous in LRLL 
31. The equivalent width is continuously changing with 
a range of a factor of 3.4, suggesting that the accretion 
rate is always varying as was found with the infrared hy- 
drogen recombination lines. Similar to the dust emission, 
the variations in the accretion rate onto the star are not 
a single isolated event, and are always occurring with 
varying amplitude. Our high-resolution Ha spectra are 
shown in Figure ITT] They show broad profiles, consistent 
with the presence of accretion. It appears as though the 
red side of the feature experiences the largest variations, 
while the blue side stays relatively constant. This may 
be due to partial obscuration of the accretion flow by the 
disk, or a change in the magnetospheric geometry. 

The outflow of material, as traced by the blue-shifted 
emission/absorption in the Hel line (Figure [TU|) . changes 
on daily to weekly timescales. In 2008 it shows large vari- 
ations from night to night, while in 2009 the line slowly 
goes from showing very strong blue-shifted absorption to 
no blue-shifted features. Blue-shifted absorption in the 
Hel line is believed to come from a combin ation of stellar 
wind and disk winds ([Fischer et al.ll2008D and variabi l- 
ity in the lines is sometimes seen ( Edwards et al.ll2006D . 
The red side of the Hel line is connected with the accre- 
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tion flow, and the strength of the red-shifted absorption 
is broadly correlated with the strength of the Pa^ and 
Br7 emission lines. When the Pa/3 hne is the strongest 
in 2005 and 2006, the red-shifted absorption is the deep- 
est, and in 2009 the absorption increases with the ac- 
cretion rate. This red-shifted absorption may result if 
we are looking down the accretion column, so the gas 
that is flowing onto the star absorbs the emission from 
the hot column of shocked gas striking the surface of 
the star. The variations in the relative strengths of the 
two sides of the line suggest that the two do not change 
contemporaneously. There are days with just red-shifted 
absorption, days with just blue-shifted absorption and 
days with both. 

There is also evidence for correlated changes in the ac- 
cretion rate and the infrared excess. In 2009 we were 
able to obtain near-infrared spectra during our Spitzer 
warm mission monitoring (vertical tick marks in Fig- 
ure [7]). We observed two spectra during the minimum 
in the light curve and three spectra after the infrared ex- 
cess had substantially increased. We find that while the 
infrared excess is at its minimum, the Pa/? line (8,9 Oct) 
is weaker than in any other epoch. As the infrared excess 
increases over the next few weeks the Pa/? line emission 
also increases. We have one observation of the Ha line 
taken on 12 Oct 2009, close to the infrared minimum, and 
the equivalent width is lower than seen in other epochs. 
This suggests that the dust at the inner edge of the disk 
and the gas flowing onto the star are closely connected 
by the physical mechanism that causes this variability, 
which has been p reviously considered a s a source of in- 
frared variability (jCarpenter et al.ll20Ql[ ). Also, the blue- 
shifted absorption in the Hel line is strongest at the in- 
frared minimum and slowly disappears as the excess gets 
stronger, suggesting that the dust and outflowing gas are 
also connected. 

4. THEORETICAL IMPLICATIONS 

Our observations of LRLL 31 target the inner disk, 
as shown by the variability timescales of a few days to 
a week and the temperature of the dust close to the 
sublimation temperature. Studies of circumstellar disks 
have gained a reasonably complete and consistent picture 
of the behavior of the accretion/protoplanetary disks at 
radii greater than 1 AU, but the zone within this radius 
presents many puzzles. Not only is it beyond our abil- 
ity to resolve well observationally, but it is the venue for 
complex processes involving intense magnetic fields, gas 
flow, accretion luminosity, and extreme heating and sub- 
limation of dust. The variability observed here can help 
to provide constraints on some of the physics occurring 
in this part of the disk. In this section, we consider pos- 
sible physical causes for this behavior such as variable 
accretion, perturbations by a companion, winds, and the 
infiuence of magnetic fields. Schematic diagrams of the 
various models considered are shown in Figure 1121 

4.1. Variable Accretion 

In a disk in hydrostatic equilibrium, the scale height of 
the dust is set partly by the lu minosity illuminating its 
surface (|Dullemond et alll2001[ ). which is a combination 
of stellar luminosity and accretion luminosity. A sub- 
stantial change in either would lead to a variation in the 
scale height of the disk (Figure [T2k ) . Our observations 



rule out any large changes in the star, but do indicate 
large variations in the accretion. Since the accretion lu- 
minosity is much smaller than the stellar luminosity, a 
much larger change in the accretion luminosity than the 
observed variation would be needed to produce the ob- 
served fluctuations in scale height. We can then rule out 
a variable illumination of the inner disk as the source of 
the variability. 

The viscous nature of the disk also means that the 
accretion rate through the disk is directly tied to its sur- 
face density (D'Alcssio et al. 1998). If the surface den- 
sity of the disk were to suddenly increase, possibly due 
to a large inflow of material from the outer disk, then 
the height of the r = 1 surface would increase (Figure 
[T2b). Outgoing stellar photons will be absorbed by the 
disk below the height of the r = 1 surface of the disk. 
A change in the height of this surface would change the 
amount of stellar light that is absorbed by the disk. This 
would lead to a change in the covering fraction of the in- 
ner disk and the shadowing of the outer disk, consistent 
with the behavior of this system. The timescale for a 
large mass of dust to travel through the inner disk is 
the viscous timescale, which is roughly 2000 years at 0.2 
AU. This is much longer than the dust continuum vari- 
ability timescales we observe. Also, if the magnetosphere 
radius where the stellar accretion flow is launched is sig- 
nificantly less than the dust sublimation radius (see be- 
low), then the long viscous timescale is also inconsistent 
with the observed correlation between dust and stellar 
accretion rate variations. 

These arguments suggest that the variations are not 
caused directly by changes in the accretion rate. It is 
possible that the relationship between the infrared excess 
and accretion is not a causal one. The change in accretion 
rate may not directly lead to a variation in the covering 
fraction of the inner disk, but instead they may both be 
effects of another physical process. 

4.2. Winds 

A highly variable wind could lead to a loss of material 
in the inner disk. This wind c ould carry both gas and 
dust out from the inner disk (jVinkovfc fc Jurkid [20071 : 
iKonigl fc Pudritd [200Q) and this rapid removal of ma- 
terial would lead to a drop in both the infrared excess 
and the accretion rate, if gas were removed just before 
it began its free fall along the magnetic field lines (Fig- 
ure 112b ). Disk wind models usually start the wind with 
a velocity on the order of the sound speed, and then 
it rapi dly accelerates to larger than the Keplerian ve- 
locity (jBlandford fc Pavnill982D . which would imply a 
slower decrease in the disk emission than is observed. 
The strong blue-shifted absorption on the Hel line during 
the minimum of the infrared excess, and the rapid fluc- 
tuations in this component are consistent with a strong 
wind. If we assume that the inner disk is just barely 
optically thick before the wind starts to remove mate- 
rial with a gas to dust mass ratio of 100, and the wind 
removes the entire inner disk in only 5 days (the length 
of the drop in the infrared light curve) then the outflow 
rate is 2 x 10~^MQ/yr. The accretion rate required to 
refill the disk over the next 30 days is 3 x 10~^MQ/yr. 
This accretion rate is a lower limit since the disk is not 
just barely optically thick but it is still much larger than 
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is observed for LRLL 31. Also, filling the inner disk by 
accretion from the outer disk should occur on a viscous 
timescale, which is much longer than the timescale for 
the increase in the infrared flux. 

The wind could be asymmetric about the plane of 
the disk, leading to a variation in the disk structure 
(Figure [T^Jl). This can arise when a stellar magnetic 
field that is not syr nmetric about the eq uatorial plane 
launches the wind ([Lovelace et al.l [2"010[) . If the mass 
flow in the wind is larger on one side of the disk than 
the other, then the asymmetric loss of angular momen- 
tum would lead to warping of the disk, which would 
the n shadow the outer disk, lead ing to infrared variabil- 
ity (|Flahertv fc Muzerolldl2010D . If the fluctuations in 
the wind were rapid enough, as is suggested by the fast 
changes in the blue-shifted side of the Hel line, then they 
could cause the observed variability. This model suggests 
that a large wind would be associated with more warping 
in the disk, and a larger infrared flux, which is inconsis- 
tent with our observations. As with the accretion flow, it 
is difficult to get a wind to explain the infrared variations 
directly, and the fluctuations seen in the blue-shifted side 
of the Hel line may be a result of the same physical pro- 
cess that changes the infrared excess and the accretion 
rate rather than the source of this behavior. 

4.3. Perturbation by a Companion 

A companion on an orbit that is misaligned with 
the disk ( Figure [T2|)__can lead to a warped disk 
(jLarwood fc FapaloizoulriQQTl ). The height of the warp 
will periodically vary as the companion passes out 
of the midplane of the disk and drags dust with it 
(|Fragner fc NelsonI [20091) . If a companion were sitting 
inside the inner disk then it could create large perturba- 
tions while also periodically disturbing the accretion rate, 
assuming it did not get close eno ugh to the inner disk to 
comp letely remove all of the dust ([Artvmowicz fc Lubowl 
Il996f ). It is also possible for a companion to be sitting 
outside of the inner disk, again assuming it isn't close 
enough to remove the inner disk. The deflcit of flux 
around 10/im relative to a typical T Tauri star indicates 
that the disk is missing material within a few AU of the 
star, i.e. it is a transition disk. Our near-infrared obser- 
vations indicate that there is still an optically thick disk 
at the dust destruction radius, but beyond this point 
there may be a gap in the disk. If so, gas and dust would 
still flow through this gap as evidenced by the continu- 
ous accretion of gas onto the star. If the gap were caused 
by a companion clearing out material near its orbit, then 
the accretion rate through the gap would be expected to 
be periodic, and material would cross the gap at a speed 
much higher than the typical velocity through a full disk 
([Artvmowicz fc Lubowlll996l ) . If the inner disk is narrow 
enough, then the shock heating as this material strikes 
the inner disk at its outer edge could lead to heating of 
the inner edge of the disk, which would change its scale 
height. 

Our data allow us to provide some constraints on the 
position and mass of a companion. A companion larger 
than IMq would be noticeable in our near-infrared spec- 
tra as anomalously strong photospheric lines in the K 
band. An equal mass binary would not produce anoma- 
lous line strengths, but would make the star appear over- 
luminous, which is not seen, although there is a large 



uncertainty in the luminosity. A smaller companion near 
the dust destruction radius can still lead to a large change 
in the structure of the dust, and would be consistent 
with the observed changes in the radial velocity. Fig- 
ure [T3| shows our constraints on the mass and location 
of a companion in the disk. The upper and lower lim- 
its on the mass come from the near-infrared spectra and 
the RV data. The limits on the location of a compan- 
ion come from the requirement that the companion can- 
not remove material from the dust destruction radius 
(taken to be 0.15 AU) or from the outer disk (^ 7 AU) 
([Artvmowicz fc Lubowl ll996). The exact location of the 
inner disk, along with the limits on the location of a com- 
panion near this dust, depends on the dust composition, 
the luminosity of the central source and the assumption 
that t he central source is a single star rather than a bi- 
nary (Na gel et aLl l2010V Moving the inner disk in or out 
will change the absolute position of the boundary, but 
not the relative positions, which are at r/a~ 0.3, 1.7 for 
the outer and inner limits respectively, where a is the 
binary separation and r is the location of the dust. The 
radial velocity data suggests that this system is a binary, 
and the presence of dust at certain locations within the 
disk limits the location of this companion. 

Even if this system is a binary, the companion may not 
be responsible for the infrared variability. A companion 
in a circular orbit in the plane of the dis k will not caus e 
substantial mid-infrared variability (Nag el et al.l l2010(). 
Making the companion misaligned with the disk could 
lead to periodic perturbations (Fragner_&_Nelson 2QD9) . 
Figure \W\ shows the periodigram based on the monitor- 
ing at 3.6 and 4.5/im. We are most sensitive to periods in 
the range 2-25 days and no signiflcant peak is seen indi- 
cating that there is no periodicity over these timescales. 
If a companion were located 0.04-0.2 AU from the cen- 
tral star and perturbed the disk on every orbit then we 
would see evidence of this in the IRAC data. This ex- 
cluded region is marked in left-hand side of Figure 1131 
This data suggests that it is unlikely that there is a com- 
panion on an inclined orbit within the inner disk, but 
does allow an object further out in the disk. Based on 
the limits discussed above, if the inner disk is at 0.15 
AU then the companion is restricted to outside of 0.4 
AU where the orbital period is ~ 100 days. The rapid 
decline in the flux at the beginning of our infrared mon- 
itoring is difficult to explain in the context of the slow 
variations expected from a companion this far from the 
star. It is possible that the inner disk is closer than 0.15 
AU, pushing in the outer limit on the position of a com- 
panion inward resulting in a timescale more comparable 
to what is observed. The companion could also be on an 
eccentric orbit, which would allow it to rapidly disturb 
the disk during its high velocity periastron motion, while 
still having a long period. 

4.4. Magnetic Fields 

A magnetic field threading a turbulent disk could 
lead to vari able structure . Using magneto-hydrodynamic 
simulations I Turner et al.l ([20101 ) find rapid variations in 
the disk due to the turbulence caused by the magneto- 
rotational instability (Figure \TT[ ) . They see sharp drops 
in the scale height of the disk by 30% in an orbital 
timescale as dust is lifted out of the disk by a buoyant 
magnetic field. There is evidence for a small change in 
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the silicate feature (jMuzeroUe et alj [20091 ) . which could 
be explained by the limiting of dust out of the midplane 
(| Turner et al.ll20Tol) . The speed of the changes in this 
model is consistent with our observations, but the small 
amount of excess in the minimum of the 3.6,4.5/im light 
curve would require that the entire inner disk have a 
scale height close to zero. Such global behavior is not 
seen in these models and the changes in the scale height 
are smaller than have been observed 

The interaction between the disk and stellar magnc- 
tosphere can be highly dynamic (Bpuvier et al. 2007b). 
According to some models, where the magnetic field 
truncates th e disk oscillates (Figure [T^ ) on timescales 
of - 25days (|Goodson fc WingleelUggaT which is similar 
to the timescale seen in the warm-mission photometry. 
The oscillations can be very large (a factor of 3 change 
in the size) due to the expansion, opening and reconnec- 
tion of the magnetosphere, although we do not see a large 
chang e in the location of the dusty disk. iRomanova et al.l 
(|2009( ) find large variations in both the accretion rate and 
the outflow rate associated with these magnetic field os- 
cillations. The timescale between bursts is ^ 60 days for 
rapidly rotating s tars, which is consistent with our in- 
frared light curve. iKulkarni fc Romanoval (|2008| ) present 
a set of three-dimensional simulations of magnetospheric 
accretion and they show that strong instabilities develop 
in the inward flow of material, breaking it up into a small 
number of tongue-shaped features extending from the 
disk, along the magnetic field lines, to the stellar surface, 
where they create hot spots (Figure fT2h). Variations in 
the number and geometry of the tongues occur on the 
dynamical timescale of the inner disk. If the tongues 
originate at or just beyond the sublimation radius and 
dust is carried along with the gas flow, they might pro- 
vide a framework for explaining our observations. They 
would carry the dust out of the disk plane, potentially 
shadow the outer disk from the star when this occurs (as 
needed to explain the longer w a velen gth variations re- 
ported in iFlahertv fc Muzerolld ()201Q) ). and would ap- 
pear as a large increase in disk scale height, all occurring 
on a timescale compatible with that of the observed vari- 
ations. 

Another possibility related to the interface of the stel- 
lar magnetic field and the disk is if the magnetic field is 
misaligned with respect to the plane of the disk (Figure 
[l2l). As material flows onto the field lines it is pulled 
out of the midplane, preferentially on one hemi s phere , 
creating a warp in the inner disk. iLai fc Zhan j poM ) 
develop a model in which the magnetic field coupling 
between the star and disk can excite waves in the disk 
when the stellar rotation axis, stellar magnetic axis and 
disk axis are misaligned. This mechanism could set in 
place a disk geometry simila r to the warps proposed by 
IFlahertv fc Muzerolld (|2010l ) to explain the variability, 
assuming either the density of the inflowing material or 
the structure of the magnetic fleld rapidly varies in order 
to change the height of the warp. 

These types of interaction require that the dust extend 
close to where the disk is truncated by the stellar magne- 

^ If the dips in the extinction that are sometimes seen in LRLL 
31 are proven to be real, the large error bars on Ay mean they 
are still marginally consistent with each other, and t he extinction 
arises from a nearly edge-on disk then the model of [Turner et al.l 
l|201CII ) would be a promising explanation for them. 



tosphere. We can estimate where material stops flowing 
in through the disk and is loaded onto the magnetic field 
lines using the corotation radius, 0.05 AU based on the 
optical photometric period. Material needs to be within 
the corotation radius before it is loaded onto the mag- 
netic field lines or else it will b e flung out in a wind ra ther 
than accreted onto the star (|Bouvier et al.ll2007b[ ). al- 
though the exact location where the magnetic field trun- 
cates the disk depends on the strength of the magnetic 
field and the accretion rate. Our observations of the in- 
frared excess find that the dust extends to the sublima- 
tion point, but translating this information into a posi- 
tion requires detailed knowledge of the density and com- 
posit ion of the dust (jlsella fc Nattal l2005t iKama et al.l 
l2009f ). Based on different assumptions about the dust 
properties and the maximum temperature of the dust 
results in a range of dust sublimation radii from 0.05 - 
0.3 AU given t he measured luminosity of L RLL31 (see 
the review by iDuUemond fc Monnieii l2010l for details 
on various calculations of the location of the dust in- 
ner rim). The smallest radii correspond to large grains 
(~ a few fim) at 1800 K, which because of their rela- 
tive high efficiency in radiating away thermal emission 
compared to small grains (O.l/im) can survive closer to 
the star before sublimating, while the largest radii cor- 
respond to small grains at 1500 K. We have no direct 
measurements of the grain properties of the inner disk, 
but the strong silicate emission feature at 10/im is consis- 
tent with the presence of small amorphous grains rather 
than the large grains needed to get the disk down to 0.05 
AU. Our measurements of the covering fraction for the 
inner disk find that it is consistently smaller than ex- 
pected for a typical T Tauri star, which could be due 
to significant increase in the grain size for the inner disk 
versus the outer disk. Substantial grain growth in the 
inner disk would move the dust edge inward, making it 
more susceptible to oscillations of the stellar magneto- 
sphere. We would expect to see variability from this pro- 
cess more often around systems with a slow rotation pe- 
riod, hence a much larger corotation radius and it would 
not preferentially occur around transition disks for which 
the type of infrared v ariability observed here is common 
(|Espaillat et al.ll2"011|) . Without detailed information on 
the exact termination of the dusty disk, it is difficult 
to determine if oscillations in the stellar magnetosphere 
cause the observed variability. 

5. CONCLUSION 

We present results from a large synoptic study of the 
transition disk variable LRLL 31. The star itself stays 
relatively constant showing variability consistent with ro- 
tation of cool spots across its surface. We find large vari- 
ations in the infrared excess and the accretion rate on 
timescales of weeks. The dust appears to stay at the dust 
destruction radius, while its scale height rapidly fluctu- 
ates by a factor of eight and the infrared flux changes by 
0.3 mag. The accretion rate varies by a factor of 5, and 
appears to be correlated with the strength of the infrared 
excess. However, the change in accretion rate is unlikely 
to be the direct cause for the change in infrared excess. 
It is also unlikely that the observed variations are due to 
the influence of a companion within 0.4 AU, based on the 
lack of periodicity in our Spitzer warm mission monitor- 
ing and the presence of an optically thick inner disk. The 
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source of the variability may be related to a companion 
outside of, but still close to, 0.4 AU or to the dynamic 
interface between the disk and the stellar magnetic field, 
assuming the dust extends close enough to the star. Our 
observations are able to limit the list of plausible physi- 
cal models, although they cannot exactly define what is 
happening to the disk. Further constraints on the radial 
velocity variations, and the exact position of the inner 
disk will help to select between these models. 

We would like to thank Paul Smith for obtaining 
the spectropolarimetric data and Joan Najita for use- 
ful discussion regarding the physical mechanism behind 
the variability. Our results are partly based on ob- 



servations obtained at the Centro Astronoomico His- 
pano Aleman (CAHA) at Calar Alto, operated jointly 
by the Max-Planck Institut fiir Astronomie and the In- 
stituto de Astrofisica de Andaluci'a (CSIC). The obser- 
vations were supported by OPTICON. OPTICON has 
received research funding from European Community's 
Sixth Framework Programme under contract number 
RII3-CT-001566. This work is based in part on obser- 
vations made with the Spitzcr Space Telescope, which 
is operated by the Jet Propulsion Laboratory, California 
Institute of Technology under a contract with NASA. 
This work was partially supported by contracts 1255094 
and 960785 from JPL/Caltech to the University of Ari- 
zona. 
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TABLE 1 
Observing Log 



Date 


MJD 


Wavelength 


Resolution 


Note 




Spex 


29 Dec 2005 


53734.2 


0.8-5;um 


~1500 




09 Nov 2006 


54049.3 


0.8-5/im 


~1500 




10 Oct 2008 


54750.5 


0.8-5/im 


~1500 




11 Oct 2008 


54751.5 


0.8-5Atm 


~800 




18 Oct 2008 


54757.4 


0.8-5Atm 


~1500 




19 Oct 2009 


54758.4 


0.8-5/im 


~1500 




08 Oct 2009 


55112.6 


0.8-2.5/xm 


~800 




09 Oct 2009 


55113.6 


0.8-2.5Atm 


~800 




31 Oct 2009 


55137.3 


0.8-5^4111 


~800 




04 Nov 2009 


55141.5 


0.8-5(Um 


~800 




08 Nov 2009 


55144.5 


0.8-5;im 


~800 




Optical Spectra 


28 Nov 2006 


54067 


4700-9000A 


~33,000 


Keck HIRES 


29 Feb 2008 


54526 


6450-6650A 


~30,000 


MMT Hectochelle 


03 Dec 2008 


54803 


4700-9000A 


~33,000 


Keck HIRES 


13 Dec 2007 


54447 


4000-7500A 


~600 


SPOL 


16 Dec 2007 


54450 


4000-7500A 


~600 


SPOL 


06 Oct 2008 


54745 


4000-7500A 


~600 


SPOL 


09 Oct 2008 


54748 


4000-7500A 


~600 


SPOL 


28 Oct 2008 


54767 


4000-7500A 


~600 


SPOL 


24 Nov 2008 


54794 


4000-7500A 


~600 


SPOL 


12 Oct 2009 


55116 


6000-9500A 


~3500 


CAFOS 


IRS 


09 Oct 2007 


54382 


5-40/im 


~600 




16 Oct 2007 


54389 


5-40/im 


~600 




24 Feb 2008 


54520 


5-40/im 


~600 




02 Mar 2008 


54527 


5-40/^m 


~600 




MIPS 


21 Feb 2004 


53056 


24;im 


photometry 




19 Sep 2004 


53267 


24/^m 


photometry 




23 Sep 2007 


53466 


24/im 


photometry 




24 Sep 2007 


54367 


24/im 


photometry 




25 Sep 2007 


54368 


24/im 


photometry 




26 Sep 2007 


54369 


24/im 


photometry 




27 Sep 2007 


54370 


24/im 


photometry 




12 Mar 2008 


54537 


24/^m 


photometry 




19 Mar 2008 


54544 


24/im 


photometry 




IRAC 


11 Feb 2004 


53046 


3.6,4.5,5.8,8.0^4111 


photometry 


Cold-mission 


08 Sep 2004 


53257 


3.6,4.5,5.8,8.0Atin 


photometry 


Cold-mission 


19 Mar 2009 


54910 


3.6,4.5,5.8,8.0Atm 


photometry 


Cold-mission 


20 Mar 2009 


54911 


3.6,4.5,5.8,8.0/im 


photometry 


Cold-mission 


21 Mar 2009 


54912 


3.6,4.5,5.8,8.0Atm 


photometry 


Cold-mission 


22 Mar 2009 


54913 


3.6,4.5,5.8,8.0/im 


photometry 


Cold-mission 


23 Mar 2009 


54914 


3.6,4.5,5.8,8.0/im 


photometry 


Cold-mission 


03 Oct-07 Nov 2009 


55107-55142 


3.6,4.5;im 


photometry 


Warm-mission 
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TABLE 2 
Photometry 



Date 


J 


H 


K 


2MASS 


12.09 


10.54 


9.69 


10 Oct 2008 


12.19 


10.71 


9.95 


11 Oct 2008 


12.23 


10.81 


10.08 


18 Oct 2008 


12.26 


10.81 


10.07 


19 Oct 2008 


12.25 


10.79 


9.99 


08 Oct 2009 


12.34 


10.92 


10.29 


09 Oct 2009 


12.22 




10.25 


31 Oct 2009 


12.26 


10.83 


10.14 


04 Nov 2009 


12.24 


10.77 


10.03 


08 Nov 2008 


12.15 


10.68 


9.92 



Note. — Typical uncertainties in the photometry, dominated by the absolute uncertainty in the flux measurements, are 0.1 mag. 



TABLE 3 

IRAC PHOTOMETRY 



Date 

19 Mar 2009 

20 Mar 2009 

21 Mar 2009 

22 Mar 2009 

23 Mar 2009 



[3.6] 



[4.5] 



[5.8] 



9.058 ± 0.002 
9.064 ± 0.002 
8.995 ± 0.002 
8.968 ± 0.002 
9.071 ± 0.002 



8.781 ± 0.003 
8.721 ± 0.003 
8.736 ± 0.003 
8.639 ± 0.003 
8.695 ± 0.002 



8.404 ± 0.006 
8.403 ± 0.006 
8.347 ± 0.006 
8.343 ± 0.006 
8.342 ± 0.006 



m 

7.997 ± 0.03 
8.033 ± 0.03 
8.002 ± 0.03 
7.941 ± 0.03 
7.981 ± 0.03 



TABLE 4 
Extinction 



Date 


Av 


Wavelength Range 


13 Dec 2007 


8.8 


Optical 


06 Oct 2008 


9.2 


Optical 


09 Oct 2008 


8.9 


Optical 


10 Oct 2008 


9.5 


NIR 


11 Oct 2008 


8.2 


NIR 


18 Oct 2008 


9.2 


NIR 


19 Oct 2008 


9.2 


NIR 


28 Oct 2008 


8.5 


Optical 


24 Nov 2008 


8.1 


Optical 


08 Oct 2009 


8.8 


NIR 


09 Oct 2009 


8.3 


NIR 


31 Oct 2009 


8.9 


NIR 


04 Nov 2009 


8.7 


NIR 


08 Nov 2009 


8.8 


NIR 



Note. — The uncertainty in the extinction is 0.5 mag. 



TABLE 5 
Stellar Luminosity 



Date Stellar Luminosity (L/Lq) 

10 Oct 2008 5.3 ± 0.8 

11 Oct 2008 3.6 ± 0.6 

18 Oct 2008 4.6 ± 0.7 

19 Oct 2008 4.6 ± 0.7 

08 Oct 2009 3.8 ± 0.6 

09 Oct 2009 3.8 ± 0.6 
31 Oct 2009 4.2 ± 0.7 
04 Nov 2008 4.1 ± 0.6 
08 Nov 2009 4.6 ±0.7 



Note. — 



Errors include the 0.1 mag uncertainty in the photometry and the 0.5 mag uncertainty in Ay- 
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TABLE 6 
Infrared Excess 



Date 


Veiling 


Dust Temperature 


Covering Fraction 


29 Dec 2005 


0.33 ±0.1 






09 Nov 2006 


0.20 ±0.1 






10 Oct 2008 


0.20 ±0.1 


1940 ± 250 


0.039 


11 Oct 2008 


0.32 ±0.1 


2100 ± 250 


0.084 


18 Oct 2008 


0.17 ±0.1 


1820 ± 200 


0.037 


19 Oct 2008 


0.16 ±0.1 


1620 ± 200 


0.031 


08 Oct 2009 


-0.05 ±0.1 




O.OIO'^ 


09 Oct 2009 


0.02 ±0.1 




0.012=* 


31 Oct 2009 


0.10 ±0.1 


1540 ± 250 


0.025 


04 Nov 2009 


0.23 ±0.1 


2020 ± 300 


0.048 


08 Nov 2009 


0.29 ±0.1 


1780 ± 250 


0.054 



Derived assuming a dust temperature of 1830 K and fitting to the IRAC photometry 



TABLE 7 

POLARIMETRY 



Date 



13 Dec 2007 
16 Dec 2007 
06 Oct 2008 
09 Oct 2008 
28 Oct 2008 
24 Nov 2008 



P(%) 



8.44 ± 0.14 
7.66 ± 0.25 
8.34 ± 0.12 
7.70 ± 0.07 
7.84 ± 0.11 
7.77 ± 0.23 



q(°) 



146.6 ± 0.5 
144.0 ± 0.9 

145.7 ± 0.4 

144.9 ± 0.3 
144.4 ± 0.4 

143.2 ±0.8 



TABLE 8 

Emission Lixe Tracers oe Accretiox Rate 



Date Pa^ EW (A) Br7 EW (A) Veiling Corrected Br7 EW (A) Veihng 



29 Dec 2005 


-3.36 


± 





08 


(0 


50) 


-3 


05 


± 





09 


(0 


84) 


-2.56 


± 


0.09 


(0 


71) 


0.33 ± 


0.1 


09 Nov 2006 


-4.29 


± 





04 


(0 


37) 


-2 


55 


± 





02 


(0 


18) 


-2.08 


± 


0.02 


(0 


15) 


0.20 ± 


0.1 


10 Oct 2008 


-2.74 


± 





07 


(0 


49) 


-1 


52 


± 





04 


(0 


20) 


-1.20 


± 


0.04 


(0 


17) 


0.20 ± 


0.1 


11 Oct 2008 


-2.69 


± 





07 


(0 


45) 


-1 


27 


± 





05 


(0 


22) 


-0.80 


± 


0.05 


(0 


16) 


0.32 ± 


0.1 


18 Oct 2008 


-1.54 


± 





04 


(0 


12) 


-0 


86 


± 





02 


(0 


07) 


-0.59 


± 


0.02 


(0 


05) 


0.17 ± 


0.1 


19 Oct 2008 


-1.78 


± 





06 


(0 


30) 


-2 


11 


± 





03 


(0 


21) 


-1.84 


± 


0.03 


(0 


19) 


0.16 ± 


0.1 


08 Oct 2009 


-1.12 


± 





12 


(0 


31) 


-1 


16 


± 





04 


(0 


14) 


-1.16 


± 


0.04 


(0 


14) 


-0.05 ± 


0.1 


09 Oct 2009 


-0.74 


± 





05 


(0 


08) 


-1 


39 


± 





03 


(0 


11) 


-1.36 


± 


0.03 


(0 


11) 


0.02 ± 


0.1 


31 Oct 2009 


-2.56 


± 





04 


(0 


22) 


-2 


04 


± 





02 


(0 


14) 


-1.87 


± 


0.02 


(0 


13) 


0.10 ± 


0.1 


04 Nov 2009 


-3.31 


± 





05 


(0 


35) 


-2 


29 


± 





03 


(0 


21) 


-1.94 


± 


0.03 


(0 


18) 


0.23 ± 


0.1 


08 Nov 2009 


-3.08 


± 





05 


(0 


32) 


-2 


30 


± 





03 


(0 


18) 


-1.87 


± 


0.03 


(0 


16) 


0.29 ± 


0.1 



Note. — 



Uncertainties in parenthesis include the uncertainties in the continuum 
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TABLE 9 

Accretion Rates 



Date Flux (10-" erg s-i cm-2) log(L/Lo) log(Lacc/L0) M {IQ-^M^yr-^) 



Pa/3 


10 Oct 2008 


16.6 ± 2.7 


-3.28 ± 0.07 


-0.59 ± 0.08 


1.62 ± 0.13 


11 Oct 2008 


11.2 ± 1.8 


-3.45 ± 0.07 


-0.78 ± 0.08 


1.05 ± 0.08 


18 Oct 2008 


8.08 ± 1.3 


-3.59 ± 0.07 


-0.94 ± 0.08 


0.73 ± 0.06 


19 Oct 2008 


9.43 ± 1.6 


-3.52 ± 0.07 


-0.86 ± 0.08 


0.87 ± 0.07 


08 Oct 2009 


4.93 ± 1.2 


-3.81 ± 0.11 


-1.19 ± 0.13 


0.41 ± 0.05 


09 Oct 2009 


3.21 ± 0.6 


-3.99 ± 0.08 


-1.40 ± 0.09 


0.25 ± 0.02 


31 Oct 2009 


12.5 ± 2.0 


-3.40 ± 0.07 


-0.73 ± 0.08 


1.18 ± 0.09 


04 Nov 2009 


15.5 ± 2.5 


-3.31 ± 0.07 


-0.62 it 0.08 


1.52 ± 0.12 


08 Nov 2009 


16.1 ± 2.6 


-3.29 ± 0.07 


-0.60 ± 0.08 


1.59 ± 0.13 


Br- 


io Oct 2008 


3.10 ± 0.38 


-4.01 ± 0.05 


-0.62 ± 0.06 


1.52 ± 0.09 


11 Oct 2008 


1.61 ± 0.29 


-4.29 ± 0.08 


-0.98 ± 0.10 


0.66 ± 0.07 


18 Oct 2008 


1.32 ± 0.18 


-4.38 ± 0.06 


-1.09 ± 0.08 


0.51 ± 0.04 


19 Oct 2008 


4.45 ± 0.48 


-3.85 ± 0.05 


-0.42 ± 0.06 


2.40 ± 0.14 


08 Oct 2009 


2.04 ± 0.26 


-4.19 ± 0.06 


-0.85 ± 0.08 


0.89 ± 0.07 


09 Oct 2009 


2.35 ± 0.26 


-4.13 ± 0.05 


-0.77 ± 0.06 


1.09 ± 0.06 


31 Oct 2009 


3.81 ± 0.40 


-3.92 ± 0.05 


-0.51 ± 0.06 


2.19 ± 0.12 


04 Nov 2009 


4.28 ± 0.46 


-3.87 ± 0.05 


-0.45 ± 0.06 


2.76 ± 0.13 


08 Nov 2009 


4.63 ± 0.50 


-3.83 ± 0.05 


-0.40 ± 0.06 


3.32 ± 0.15 



Note. — Uncertainties include the uncertainty in the spectra and the J and K photometry, but do not include the uncertainty in the 
conversion from line flux to a<;cretion luminosity. The Br7 line fluxes are derived with veiling of the line taJcen into a<;count. 



TABLE 10 
Ha EW 



Date Ha EW(A) 

28 Nov 2006 11.5 

13 Dec 2007 4.86 

29 Feb 2008 12.0 
06 Oct 2008 12.06 
09 Oct 2008 10.32 
28 Oct 2008 10.77 
24 Nov 2008 13.56 
03 Dec 2008 8.12 
12 Oct 2009 4.02 
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TABLE 11 






IRAC Warm Mission Photometry 




MJD-55000 


[3.6] 


[4.5] 


106.91 


9.314 


9.098 


107.30 


9.373 


9.109 


107.57 


9.363 


9.127 


107.81 


9.437 


9.134 


108.13 


9.439 


9.174 


108.38 


9.469 


9.214 


108.58 


9.476 


9.200 


108.81 


9.482 


9.236 


109.58 


9.498 


9.271 


110.50 


9.514 


9.306 


111.07 


9.609 


9.372 


112.43 


9.642 


9.438 


113.07 


9.639 


9.413 


114.73 


9.559 


9.341 


115.37 


9.538 


9.302 


115.97 


9.552 


9.272 


116.32 


9.495 


9.277 


116.62 


9.529 


9.303 


116.72 


9.527 


9.313 


117.13 


9.535 


9.288 


117.21 


9.533 


9.299 


117.46 


9.551 


9.319 


117.69 


9.565 


9.311 


118.60 


9.568 


9.321 


119.25 


9.565 


9.325 


119.93 


9.506 


9.289 


121.21 


9.510 


9.268 


123.07 


9.491 


9.206 


124.96 


9.488 


9.282 


127.68 


9.472 


9.244 


129.81 


9.475 


9.212 


131.27 


9.454 


9.171 


133.38 


9.363 


9.097 


135.61 


9.432 


9.199 


137.52 


9.333 


9.075 


139.07 


9.349 


9.058 


141.52 


9.221 


8.937 


142.13 


9.245 


8.922 



Note. — Uncertainties are 0.017 and 0.014 mag at [3.6] and [4.5] respectively based on the rms fluctuations of non-variable stars in the 
field 
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Fig. 1. — Observed 0.8-5/im spectra of LRLL 31 (black). Also included are simultaneous JHK and 3.6,4.5/im photometry (blue), when 
available. The spectra are compared to the standard with our best fit reddening (red spectra). All spectra have been scaled to match the 
near-infrared photometry, except the 2005 and 2006 spectra for which we do not have contemporaneous photometry and have not been 
flux calibrated. For the photometric points the 'error bars' in the x-axis show the size of the passband. 
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Fig. 2. — Fig [T] continued. Curved shape of the H and K band near 1.8/im in 2006, which is not seen in the other spectra, is a sign of 
imperfect telluric correction. 
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Fig. 3. — Optical spectra of LRLL 31 normalized at 6OOOA. Red dashed line shows a reddened kurucz model that best matches the shape 
of the spectra. 
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Fig. 4. — Continuum normalized flux as a function of wavelength for the 16 bands used to measure the veiling for 08 Nov 2009. In each 
panel both the observed spectra (black line), along with the best fit veiled standard (red dashed line) are shown. The veiling increases 
from in J band to a few tenths in K band. There can be a large spread in the veiling from band to band depending on the strength of 
the lines in the band and the continuum fit. 
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Fig. 5. — The infrared excess. Each plot show the difference between LRLL 31 and the standard, normalized to the photospheric flux at 
2.15/im. On each plot the best fit is shown with a red dashed line and 3.6,4.5/im photometry is included when available. Each spectra has 
been smoothed by a median filter 0.01/im wide in order to reduce the noise in the continuum. The parts of the spectra marked in blue are 
strongly affected by the telluric correction. Note that the scale changes for each plot. The strength of the excess shows large variations 
from night to night, although the temperature stays roughly constant. 
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Fig. 6. — The spectra for each night, along with a standard and the blackbody fit the the excess. Dark blue line is the standard, red 
dashed line is the standard plus the best fit blackbody. Blue squares are the veiling measurements, and blue triangles are the 3.5,4.5^im 
photometry, where available. Spectra have been normalized by the photosphere at 2.15^m and smoothed by a median filter O.Olfim wide. 
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Fig. 7. — 3.6 and 4.5/im light curves for LRLL 3f . The data show a rapid drop in flux followed by a slow increase over the course of a 
few weeks. Uncertainties in the photometry are similar to the size of the points. Vertical lines mark the times when we obtained spectra. 
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Fig. 8. — Spectral energy distribution of LRLL 31 showing our early IRS data (rereduced), the IRAC-cryo data, the JHK photometry 
and the range of magnitudes seen in the IRAC warm mission data (error bars at 3.6,4.5/xm). A kurucz model photosphere, reddened by 
Av=8.8 and scaled to the average J band magnitude has been included for comparison. Vertical lines near the near-infrared photometry 
demonstrate the uncertainties in this data. Uncertainties in the mid-infrared photometry are similar to the size of the data points 
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Fig. 9.— Line measurements from 2005,2006, 2008 (left two panels) and 2009 (right two panels). Both the Pa/3 Une (A = 1.282/im) 
and Br7 (A = 2.16/im) are shown for each day. The solid line is the observed spectrum, while the dashed line has had the photospheric 
absorption subtracted from the line. At the top of each panel is the G6 WTTS standard. There are large variations in the emission lines, 
mainly on timescales of weeks. 
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Fig. 10. — Hel (A = 1.083/im) and Pa/3 (A = 1.0941/im) lines in LRLL 31. The Hel line shows large variations in both the blue and the 
red side in as little as one day. 
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Fig. 11. — High resolution Ha measurements, normalized to the continuum. The line shows continuous variations in its shape and 
strength. 
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Fig. 12. — Schematic diagrams of the varioijs models considered for the variability of LRLL 31. 



25 




0.01 0.10 1.00 5 10 15 20 25 30 35 

Seporotion (AU) Doys 



Fig. 13. — (Left) Limits on tlic mass of a companion. Hashed areas siiow regions of mass and position of a companion that can be 
ruled out based on our data. The upper limit is from the near-infrared spectra and the lower limit is from our RV data (dashed line 
assumes i=38°). The limits on position come from the dynamical clearing of the disk by a companion. Disk material exists at the dust 
sublimation radius (~0.15 AU, dotted line) and in the outer disk (~7 AU), and we can rule out a companion close to these locations. The 
lack of periodicity in the 3.6,4.5/im photometry rules out a companion within the grey area perturbing material on every orbit. (Right) 
Lomb-Scargle Normalized Periodigram based on the 3.6 and 4.5/im photometry. The dashed lined indicates a false alarm probability of 
0.1%. There are no peaks above this line, indicating that there is no significant periodic signal in the data. 



